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ABSTRACT 




/*. I 



A computer algorithm was developed to determine if an acoustic transmitter can 
be localized based on estimates of local angles of arrival of acoustic signals incident 
upon a receive planar sonar array, knowledge of the deterministic effects of the ocean 
on sound propagation, and local sound-speed profiles of the ocean. The algorithm was 
designed to determine azimuthal and elevation/depression angles to the transmitter as 
well as computing the depth, range, cross range, and line-of-sight range separations 
between the transmitter and the receive array. The algorithm utilizes ray acoustics and 
model-based phase weights to determine the transmitter's location relative to the 
receive array's position. As written, the algorithm is capable of solving localization 
problems in which the transmitter and receiver are in the same gradient of the local 
sound-speed profile, provided that the range from transmitter to receiver is not so great 
that the acoustic signal passes through a turning point prior to reaching the receive 
array. The results indicate that the method proposed is viable for the class of problems 
for which it was designed, and accuracies on the order of 0. 1 meters are obtained for 
line-of-sight ranges on the order of several kilometers. The angles calculated by the 
algorithm are all accurate to within 0.005 degrees. 
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I. INTRODUCTION 



This thesis constitutes one part of a long range project to develop new sonar 
signal processing algorithms capable of rapidly solving sonar localization problems. At 
present, the solution of the sonar fire control problem can require a considerably 
longer time than that required for most other types of fire control problems. The long 
time required to achieve a solution can cause a significant degradation in a ship's 
ability to avoid counterdetection, due to continuously decreasing range to the target 
during problem solution. A sonar system capable of rapid target localization without 
requiring own ship's maneuvers would greatly enhance the capabilities of our ships, 
and allow for weapon firings at longer ranges. 

The research question investigated in this thesis is whether or not it is possible 
to develop an algorithm which utilizes estimates of the local angles of arrival of 
acoustic signals incident upon a planar sonar array, knowledge of the deterministic 
effects of the ocean medium on sound propagation, and local sound-speed profiles of 
the ocean, to locate an acoustic transmitter, both in azimuthal angle and 
elevation depression angle. In addition the model-based localization algorithm 
(hereafter referred to as the 'localization algorithm') was designed to provide the range, 
depth, cross range, and line-of-sight range between the acoustic transmitter and the 
receive array. 

Ray acoustics provides methods of determining ranges and propagation angles 
for transmission of acoustic signals in inhomogeneous media [Ref. 1: sect. 6.2]. The 
deterministic effects of the inhomogeneous ocean medium on acoustic signals are well 
known. From a transmitter in a known position, it is possible to develop ray traces 
that illustrate the propagation of acoustic signals through the ocean medium. The 
intent here is to use this knowledge of sound propagation to find the transmitter's 
location based on the estimated angles of arrival at a receive array. The estimates of 
the local angles of arrival are obtained from a frequency domain adaptive beamforming 
algorithm developed by Ziomek and Chan [Ref 2 \. This algorithm performs frequency 
domain adaptive beamforming for planar sonar arrays using a modified complex LMS 
adaptive algorithm. The algorithm generates estimates of the local angles of arrival, 
namely, the azimuthal and elevation depression angles, of incoming acoustic signals. 
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However, in a real ocean environment, these local angles of arrival do not reflect the 
true line-of-sight angles to the target. 

The localization algorithm uses the angle-of-arrival estimates, plus typical 
sound-speed profiles that are normally available to ships. In addition, it was found 
that one more piece of information is required to localize the target. This information 
is a model-based phase weight which is part of a model-based signal processing 
algorithm developed by Ziomek and Blount [Ref. 3]. These phase weights are used to 
"correct for deterministic, ocean medium, phase effects due to ray bending as a signal 
propagates in the inhomogeneous ocean medium whose index of refraction (sound- 
speed profile) is a function of depth." [Ref. 3] The phase weights were originally 
developed as part of an underwater acoustic communication problem in which receiver 
and transmitter locations were known. The form of the phase weights used will be 
presented in Chapter II. 

For the problem investigated in this thesis, transmitter location is unknown a 
priori and, therefore, the model-based phase weights cannot be determined in exactly 
the same manner as was done in the algorithm developed by Ziomek and Blount 
[Ref. 3]. The usefulness of the localization algorithm developed in this thesis is based 
on the availability of the model-based phase weights. The research done here is a 
feasibility study of the ability to localize an acoustic transmitter if the phase weights 
were available. The development of an algorithm to generate the model-based phase 
weights was not the subject of this research. 

The localization algorithm is limited to solving a particular class of problems. 
The localization algorithm is designed to accommodate vertical variations in sound- 
speed profile or, sound-speed profiles that are functions of depth only. Horizontal or 
range variations in sound-speed profile were not examined in this initial study because 
they constitute only a relatively small portion of ocean areas. Additionally, the 
transmitter and receiver are assumed to both be within the same sound-speed gradient. 
Finally, all case studies were conducted based on the assumption that the receiver was 
in close enough proximity to the transmitter so that the acoustic signal had not passed 
through a turning point prior to reaching the receiver. A turning point is defined as 
the point along a ray path at which the angle of propagation is 90 degrees with respect 
to the positive Y, or depth, axis. These three restrictions were necessary to limit the 
scope of the initial study to a size that would allow for a complete verification of the 
localization technique proposed, in the time allotted for the study. 
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Chapter II describes the theory used to develop the localization algorithm. An 
overview of the problem and its geometry is presented, and then the computations 
leading to the algorithm are discussed. Finally, the limitations of the algorithm are 
presented. 

Chapter III consists of the computer simulation results and an explanation of the 
implementation of the theory in a computer algorithm. The output from the 
localization algorithm is compared to the known geometry, and a comparison of 
double precision versus single precision results is included. Additionally, the program 
is investigated to determine if errors develop as a function of the transmission angle 
and or depth separation. As will be shown in Chapter II, the roots of a fourth-order 
polynomial must be determined to find the angle of transmission at the source. The 
roots for the fourth-order polynomial are found through use of an International 
Mathematical Subroutine Library (IMSL) subroutine and are verified by comparison 
with graphs of the function. These graphs also assist in determining the correct root to 
use during problem solution. 

In Chapter IV, conclusions and recommendations are presented. 
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II. THEORY 



A. PROBLEM OVERVIEW AND GEOMETRY 

Traditionally, the localization of acoustic transmitters by ships has been carried 
out by obtaining many lines of bearing to the transmitter, and comparing these with 
own ship's motion to develop a geographic picture of the transmitter's motion. This 
method is time consuming and usually very lacking in terms of accuracy. Due to the 
nature of the deterministic effects of the ocean medium, a great deal of information is 
contained in the angles at which acoustic energy arrives at the receiver. Extraction of 
this information from the local angles of arrival, while not a simple task in of itself, 
would greatly simplify the problem of target localization. 

As a first step in exploiting the information contained within the local angles of 
arrival, a geometry must be assumed for the problem. Figure 2.1 illustrates the general 
three-dimensional geometry used in the development of the method of target 
localization presented here. 

From Figure 2.1 the following definitions are apparent: 

rectangular coordinates of the transmitter in meters, 
rectangular coordinates of the center of the receive planar 



• x 0’ y 0 > z o 

• x R . y R , zr 

• AX, AY, AZ 



• AR 



• HDLTR 



array in meters. 

cross range, depth, and Z coordinate separations, respectively, 
in meters between the transmitter and the receive array, 
where: 

AX = x R - x Q 
AY = y R - y Q 
AZ Zr - z.q 

polar radial distance in meters from the transmitter to the 
receive array. 

Note: AR 2 = AX 2 + AZ 2 



polar radial distance in meters that a ray would travel in a 
homogeneous medium (constant sound-speed profile) between 
depths y Q and Vr based on an angle of transmission of p(y Q ). 

• HDLTX, HDLTZ distances in the X and Z directions, respectively, that a ray 

would travel in a homogeneous medium between depths v Q 
and Vj^ based on an angle of transmission of p(y Q ). 

• HRLOS line-of-sight range that a ray would travel in a homogeneous 

medium between depths y Q and \r based on an angle of 
transmission ofp(y Q ). 
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Figure 2.1 System Geometry. 
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(DEPTH) 



Note: HRLOS 2 = HDLTX 2 + AY 2 + HD LIZ 2 . 

• RLOS line-of-sight range between the transmitter and the center of 

the receive array. 

Note: RLOS 2 = AX 2 + AY 2 + AZ 2 

• P(Vq) initial angle of propagation (angle of transmission), measured 

with respect to the positive Y axis, of the acoustic signal at 

source depth y 0 meters. 

• P(vj^) angle of arrival of incident plane wave field at depth 

meters. 

• PLOS the line-of-sight angle, as measured from the positive Y axis, 

between the transmitter and the receive array. 

Note that in Figure 2.1 the positive Y axis is defined in the direction of increasing 
depth, or in the downward direction. The coordinate system shown in Figure 2.1 is 
applicable for any relative positioning of the transmitter and receive array, even if AX, 
AY, and/or AZ are negative. Thus, the algorithm will work for any direction of arrival 
of the incident acoustic plane-wave field. 

The receive array is assumed to possess knowledge of its own depth. In addition, 
the receive array will have available estimates of arrival direction cosines associated 
with the local angles of arrival. These estimates are computed by the frequency 
domain adaptive beamforming algorithm. From these known quantities and 
information about the local sound-speed profile, the transmitter's location with respect 
to the receiver shall be determined. 

B. TRANSMITTER LOCALIZATION THEORY 

Energy, whether it is acoustic or electromagnetic, will refract as it passes from a 
medium with index of refraction n t into a medium with index of refraction n 2 , provided 
that n t * n 2 . In this study, the ocean volume is characterized by a one-dimensional 
index of refraction (sound-speed profile) that is a function of depth. Snell's law is 
given by [Ref. 1: p. 218], 

sin P(y) = sin P(y 0 ) 

<-’(}’) c(y Q ) 

where c(y) is the speed of sound in meters per second at a depth y. From Snell's law a 
ray parameter may be defined as 
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( 2 . 2 ) 



sin P(y 0 ) = sin p(y R ) = sin P(y TP ) = 1_ 

c(y 0 ) c(y R ) c(y Tp ) c(y-j-p) 



where: 

• b is the ray parameter. 

• yjp is the depth of a turning point. A turning point is defined as the point 

along a ray path at which the angle of propagation, P(}yp), is equal to 
90 degrees. 

At this point p(y R ) is known, since the direction cosine 

v(y R ) = cos P(y R ) (2.3) 

is calculated by the frequency domain adaptive beamforming algorithm. The speed of 
sound at depth y R , denoted c(yn), is normally known aboard ship as a result of 
measurements made by onboard sonar systems. 

It is assumed that the sound-speed profile is a linear function of depth with 
constant gradient. In most areas of the ocean this is a good approximation if both the 
transmitter and the receive array are in the same portion of the sound-speed profile. A 
typical sound-speed profile is shown in Figure 2.2. The parameter g is the constant 
gradient of the sound-speed profile in seconds' 1 . From the surface to about 100 meters 
a positive gradient is typically observed with a gradient g ~ +0.016 sec' 1 
[Ref. 4: p. 30], [Ref. 5: p. 401]. Below 100 meters a negative gradient is present, and in 
this example g ~ -0.02956 sec' 1 . Finally, at depths between 700 to 1500 meters 
[Ref. 4: p. 32] the gradient reverts to a positive value of g ~ +0.017 sec' 1 
[Ref. 5: p. 401]. The value of g in the negative portion of the gradient was computed 
by assuming the speed of sound to be 1500 meters per second at the ocean surface and 
1475 meters per second at a depth of 1000 meters [Ref. 6: p. 3). A depth of 1000 
meters was chosen as the starting point of the second positive gradient. The negative 
gradient was then calculated to fit between the positive gradients. Based on the 
assumption that both the transmitter and the receive array are in the same gradient of 
the sound-speed profile, the speed of sound at depth y can be found from 



c(y) = c(v 0 ) + g(y - y 0 ) . 



(2.4) 
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0 



Speed of Sound 




Figure 2.2 Typical Sound-Speed Profile. 

The radius of curvature that describes the arc of the circle followed by an 
acoustic field propagating through this medium is then (Ref. 1: p. 237] 



R c 

All the terms on the far righthand side of Equation 2.5 are known. Equation 2.4 may 
be rewritten as expressed by Ziomck (Ref. 1: p. 238] 



L '( > n) 



c( >'r) 



|g sin P(y Q )| |g sin p(y R )| 



(2.5) 



y 



+ 



c(y) - c(y 0 ) 
g 



(2.6) 



Therefore, 
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(2.7) 



1 1 

AY = ?R * ?o = — c(y R } ' T c(y o ) 



and, from equations 2.1 and 2.2, 



c (v R ) = 



sin P(y R ) 



( 2 . 8 ) 



and 



c (>n) = 



sin P(y n ) 



(2.9) 



Combining equations 2.8 and 2.9 with equation 2.7 it is readily observed that, 



= >'R - >’o = ^ sin P( y R) - ^ sin P(y 0 ) 



(2.10) 



or, 



AY = y R - y Q = a sin p(y R ) - a sin p(y Q ) 



(2.11) 



where 




The only unknowns now in equation 2.1 1 are P(y 0 ) and AY. Also note that 

R c = |a| = radius of curvature. (2.13) 

The radial distance AR shown in Figure 2.1 can be found by utilizing the 
following equation [Ref. 1: p. 238]: 
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(2.14) 



z = z 0 + 1 -~wt : t cos P(y 0 ) ■ cos P(y)J 

g sinp (y 0 ) u 

which is the Z coordinate of a ray propagating in the YZ plane. In this thesis a more 
general class of problem is assumed so that the coordinate axes can remain fixed 
relative to the platform on which the planar array is mounted. Therefore, in a three- 
dimensional system, z and z Q are replaced by the polar coordinates r and r Q to give 



r = r 0 + 



c (>q) 

g sinp (y n ) 



[cos p(y 0 ) - cos p(y>] , 



(2.15) 



and, as a result, 



ar = r R * r o = ~ g i cos P(y 0 > ' cos P(yR)] 



(2.16) 



or 



AR = r R ■ r o = a cos P(y 0 ) * a cos P(yR) • (2.n) 

The only unknowns in equation 2.17 are P(y 0 ) and AR. Also, note in Figure 2.1 that if 
AX = 0, (2.18) 



then 

AR = AZ . (2.19) 

At this point ray acoustics cannot provide any further information to develop a 
solution to the problem. However, a model-based phase weight for a planar sonar 
array, similar to that shown in Figure 2.3 , can be used to localize the transmitter. As 
derived by Ziomek and Blount [Ref. 3[ 

O n (Q = -2;tf Y nd Y + 0 MD (f,n) n = -(N-l)/2,...,0,...,(N-l)/2 (2.20) 
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Figure 2.3 Receive Planar Array Geometry. 
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where 



C _ - v B f 

Y c ^r) ’ 


(2.21) 


V B = cos P(y T ) , 


(2.22) 



and: 



•«vo 



• f 



• yj 

• dy 



is the phase weight in the Y direction associated with element (m,n) 
in the receive array. 

is the frequency of the transmitted electrical signal. 

is the model-based phase weight which is related to the deterministic 
angle modulation performed by the ocean medium on the transmitted 
electrical signal as a function of depth [Ref. 3]. 

is the depth of the transmit array. 

is the interelement spacing in the Y direction associated with the 
receive array. 



Equation 2.20 describes the phase weights in the Y direction that a planar sonar 
array using the three-dimensional FFT beamformer presented by Ziomek and Blount 
[Ref. 3] would use to receive an acoustic signal transmitted from a depth yy and 

received at a depth y R . This equation can be seen to consist of two parts. The first 

/ 

portion is the term -27tf yndy which is the phase weight used in traditional beam 
steering. The second part, <Pyj£)(f,n), is further described by Ziomek and Blount 
[Ref. 3] as 



<P MD (f,n) = - [k(y T )/2v B ]{[c(y T )/g][n D (y R + ndy) - 1] + \YJ (2.23) 

where the wave number in radians per meter as a function of depth yy is 

k(yy) = 2*G'c(yy) , (2.24) 



f= f c + 



kf„ 



k - -K,...,0,..-.,K , 



(2.25) 
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(2.26) 



n d(>'r + nd y) = 



c(y T ) 



c (y T ) + gAY n ’ 



and, 



AY n - >'R ‘ >'T + nd Y - 



(2.27) 



where: 



• f is the carrier frequency in hertz, 

• f Q is the fundamental frequency in hertz of the finite Fourier series representation 

of the complex envelope of the transmitted electrical signal, and 

• K is the highest harmonic used in the finite Fourier series. 

The term n'j^ defines an index of refraction which is corrected for the distance 
that the (m,n) element in the receive array is offset in the Y direction from the center 
of the array. This compensation is provided by AY n , which computes the depth 
separation between the center of the transmit array and the element (m,n). 

When using these model-based phase weights it should be noted that y-j- is 
equivalent to y Q of Figure 2.1. Additionally, 



v 



B = cos P(y T ) = cos P(y 0 ) . 



(2.28) 



Dividing equation 2.24 by 2 V 3 yields 



k(y y ) _ 27tf 

2 v b c(yqv 



c(y T ) 2 \'b c(yj )\’3 



(2.29) 



The term n'j)(y^ + ndy) - 1 may also be rewritten as 



n d(Yr + nd Y^ ' 




c (>‘j) 



(2.30) 



and, as a result 
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n 'o(yR + nd y) ■ 1 



c(yj) ~ c (>’t^ ~ gAY n 

c (y T ) + gAY n 



(2.31) 



n D (y R + nd y) • 1 



-§ AY n 

c (y T ) + gAY n 



(2.32) 



Therefore, substituting equations 2.32 and 2.29 into equation 2.23 yields 



°MD( f » n ) 



-nf [-c(y T )AY n + c(y T )AY n + gAY n 2 ] 
c(y T )v B [c(y T ) + gAY n ] 






-*f gAY n 2 

c(y T )v B [c(y T ) + gAYJ ' 



(2.33) 



(2.34) 



Expanding the denominator of the second term on the right side of equation 2.34 
results in 



c (y T ) + gAY n = c(y T ) + g(y R - y T + nd y ) • 



(2.35) 



From equation 2.4 it can be seen that 



c(y T ) + g(y R - y T + nd y ) = c(y R + nd y ) . 



(2.36) 



Therefore, 



c(y T ) + §AY n - c(y R + nd y ) 



(2.37) 



Substituting equation 2.37 and equation 2.22 into equation 2.34 gives 



-ttfgAY 2 

°MD( f ’ n ) - z- 



c(y T )c(y R + nd Y ) cos P(y T ) 



(2.38) 



From equation 2.2, with y y = y Q 
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(2.39) 



c(y T ) 



sin p(y T ) 
b 



and, as a result, 



OMD( r ’ n ) 



-7ilbgAY n 2 

c(y R +nd Y ) sin p(y T ) cos p(y T ) ' 



From equation 2.12, 



a = 



1 

bg 



or 



1 




Using equation 2.41 in equation 2.40 yields 

-JtfAY 2 

Oun(f,n) — 

• V1U ac (YR + nd y) si n P(Yj) cos P(Yj) 



(2.40) 



(2.12) 



(2.41) 



(2.42) 



where 

AY n = y R - yj + ndy = AY + ndy . (2.43) 

If the center element of the receive array is chosen as the element at which the 
phase weight O^j-^f.n) is calculated, then n = 0, and 

AY 0 = AY = y R - y T . (2.44) 

Therefore, at n = 0 
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(2.45) 



-JtfAY 2 

°MD( f ’°) ac(y R ) sin p(y T ) cos p(y T ) ' 

Now squaring both sides of equation 2.10 will result in 

AY 2 = a 2 sin 2 P(y 0 ) - 2a 2 sin p(y R ) sin p(y 0 ) + a 2 sin 2 p(y R ) (2.46) 

which can be used in equation 2.45 to replace AY 2 . Now let 

x = sin P(y 0 ) (2.47) 

and 

y = cos P(y 0 ) . (2.48) 

The x and y defined in equations 2.47 and 2.48 are not the x and y coordinates 
related to Figure 2.1. Rather, this x and y are merely dummy variables to be used in 
the solution of equation 2.45. Due to the definitions of equations 2.47 and 2.48 a 
relationship between the variables x and y is apparent, that is, 

x 2 + y 2 = 1 (2.49) 

and, therefore, 

y = ± (1 - x 2 ) 1 / 2 . (2.50) 

Next, replace yj with y Q in equation 2.45, substitute equation 2.46 into equation 
2.45, and multiply equation 2.45 by ac(v R )xy. This results in 

ac(y R )0> MD (f,0)xy = -rtf [a 2 x 2 - 2a 2 sin p(y R )x + a 2 sin 2 p(y R )] . (2.51) 

Divide both sides of equation 2.51 by nfa 2 to get 
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(2.52) 



= -x 2 + 2 sin P(y R )x - sin 2 p(y R ) . 
Rewriting equation 2.52 yields 

x 2 + <J> MD (f,0)xy - 2 sin p(y R )x + sin : p(y R ) = 0 

or 

Ax 2 + Bxy + Cx + D = 0 
where 

A = 1.0 , 



B = 



c( ?r> 

rcfa 






C = -2 sin P(y R ) , 
and 

D = sin 2 p(y R ) . 

Substituting equation 2.50 into equation 2.54 yields 
Ax 2 ± Bx(l - x 2 ) 1,/2 + Cx + D = 0 



or, 



Ax 2 + Cx + D = ± Bx(l - x 2 ) l/2 . 



(2.53) 

(2.54) 

(2.55) 

(2.56) 

(2.57) 

(2.58) 

(2.59) 

(2.60) 
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Squaring both sides of equation 2.60 gives 



(Ax 2 + Cx +D) 2 = B 2 x 2 (1 - x 2 ) = B 2 x 2 - B 4 x 4 . 



(2.61) 



Expanding equation 2.61 yields 

A 2 x 4 + 2ACx 3 + (2AD + C 2 )x 2 + 2CDx + D 2 = B 2 x 2 - B 2 x 4 (2.62) 



or 



(A 2 + B 2 )x 4 + 2ACx 3 + (2AD - B 2 +C 2 )x 2 + 2CDx + D 2 = 0 . (2.63) 



To find the unknown, x, the roots of equation 2.63 must be computed. These 
roots will also be the roots of equation 2.59. In the computer algorithm written to 
implement this theory, the value 



F(x,y) = Ax 2 ± Bx(l - x 2 ) 1/ 2 + Cx + D 



(2.64) 



was also calculated to verify the validity of the roots found for equation 2.63. 

Recalling that equation 2.47 defined 

x = sin P(y Q ) (2.47) 

and equation 2.48 defined 

y = cos P(y 0 ) , (2.48) 

we see that once x and y are known they may be substituted into equations 2.11 and 

2.17 to solve for AY and AR (since P(y^), the radius of curvature (a), the receive array 
depth, and the local sound-speed profile are all known). 

At this point AY, AR and P(y 0 ) are known. The next values to be found are AX, 
AZ, RLOS, and pLOS. Using the definitions of the direction cosines as presented by 
Ziomek [Ref. 1: p.226] 

v = cos P(y), (2.65) 
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u = cos a(y), 



( 2 . 66 ) 



and 



w = cos y(y) 



(2.67) 



where: 

• a(y) is the angle at a depth y measured from the positive X axis to the 

vector of interest. 

• y(y) is the angle at a depth y measured from the positive Z axis to the 

vector of interest. 

Referring to Figure 2.1, the direction cosine v(y) at the transmitter depth can be 
written as 



v(y 0 ) = cos P(y Q ) = 



AY 

HRLOS 



( 2 . 68 ) 



and, as a result, 



HRLOS = 



AY 
v (> 0 ) ’ 



(2.69) 



Also from Figure 2.1 it can be observed that 

HRLOS 2 = HDLTR 2 + AY 2 . (2.70) 



In ray acoustics, as presented by Ziomek [Ref. 1: p.223], the propagation vector 
is defined as 

k = k^x + kyy +k 2 z (2.71) 



where 



k x - k Q u , 



(2.72) 
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(2.73) 



ky = k 0 v , 



k z = k 0 w , 
and 



k o 



2nf 

c (>n) 



Therefore, at the transmitter, 



27tf 

kxT ' «£> u(yT) ' 



and, at the receive array, 



2;tf 

kxR = u(yR) 



Additionally, for an inhomogeneous medium which has 
function of depth only, it is known that [Ref. 1: p.223] 

ky R = kyy = constant . 

Therefore, from equations 2.76 and 2.77, 



2ttf 

C(Vr) 



u(y R ) = 



2xcf 

c (>'x) 



u(y T ) 



so that 



u(y R ) = 



c ^r) 

c(y T ) 



u(y T ) 



(2.74) 

(2.75) 

(2.76) 

(2.77) 

a sound-speed profile that is a 

(2.78) 

(2.79) 

(2.80) 
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or, 



u (}'t) ~ 



c (yy) 

c (y R ) 



u(y R ) • 



(2.81) 



In equation 2.81 u(y^) is supplied by the beamformer, c(y^) is known by own 
ship, and since equation 2.11 has been solved for AY, it is possible to use equation 2.4 
to calculate c(y-p). Therefore, u(\j) becomes a known quantity. An alternate method 
of determining c(y-p) would be by the use of Snell's law, or equation 2.1, since P(y R ), 
P(vy) and c(yj^) are all known. 

Similarly [Ref. 1: p. 233], 



k ZR “ k ZT 
and, as a result, 



(2.82) 



w(y T ) = - - * w(vn) 

T c(y R ) R ' 



(2.83) 



Referring to Figure 2.1 and utilizing equation 2.81 it can be seen that 



u(y 0 ) = cos «(y Q ) = 




u(y R ) 



HDLTX 
HRLOS ‘ 



(2.84) 



Therefore, since u(y Q ) is known from substituting y Q for \j in equation 2.81, the value 
of HDLTX is given by 



HDLTX = u(y 0 )HRLOS . 



(2.85) 



Now that u(y Q ) and v(y 0 ) are known from equation 2.84 and equation 2.68, it is 
possible to find w(y Q ) by use of the fact that [Ref. 1: p. 224] 



w~(y ) = 1 - u 2 (v ) - v*(y n ) 



( 2 . 86 ) 
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Again utilizing Figure 2.1 and the fact that 



vv(y 0 ) = cos y(y 0 ) , 



we see that 



HDLTZ 

w(y n ) = . 

0 HRLOS 



(2.87) 



( 2 . 88 ) 



Therefore, 

HDLTZ = w(y 0 ) HRLOS . (2.89) 



Figure 2.4 shows the geometry of Figure 2.1 as seen by looking down into the 
XZ plane from above the transmitter's depth. From Figure 2.4 the relationships 
between AZ, AX, and AR may be derived. 

The angle 8 in Figure 2.4 can be found from 



tan 6 = 



HDLTX 
HDLTZ ' 



(2.90) 



Therefore, 

8 = tan'^HDLTX/HDLTZ) . (2.91) 

Substituting equations 2.85 and 2.89 into equation 2.91 results in 

8 = tan' 1 ([u(y 0 )HRLOS]/[w(y 0 )HRLOS]} (2.92) 

so that 

8 = tan*’[u(y 0 )/w(y 0 )] . (2.93) 
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HDLTZ 



A Z 
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► 




Ax 



Figure 2.4 Topview of Geometry. 

For an inhomogeneous medium with a sound-speed profile that is a function of 
depth only, it can be shown that [Ref. 1: p. 232J 



«(y) 

= constant . 

w(y) 



(2.94) 



Therefore, 

6 = tan' l [u(y R )/w(y R )l (2.95) 

where u(y^) and w(y^) are available from the frequency domain adaptive bcamformer. 
From Figure 2.4, AZ and AX arc given by 

AZ = AR cos 5 (2.96) 



and 
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AX = AR cos 8 . 



(2.97) 



Referring once again to Figure 2.1, RLOS is given by 
RLOS = (AX 2 + AY 2 + AZ 2 ) 1 '' 2 
or, since 

AR 2 = AX 2 + AZ 2 , 



(2.98) 



(2.99) 



RLOS = (AR 2 + AY 2 ) 1 ' 2 . (2.100) 

Finally, PLOS can be determined by using 

pLOS = cos _1 (AY/RLOS) . (2.101) 

The equations presented in this section comprise the theory used to develop the 
model-based localization algorithm. By the use of ray acoustics and the assumption 
that the model-based phase weight is known, a closed form solution is possible for the 
localization problem. Obviously, the solution's accuracy depends on a ship's ability to 
correctly measure the sound-speed profile and the effects of any other local sonar 
conditions, such as shallow depths and the presence of biologies. However, in the open 
ocean, when the transmitter and receiver are located in the same gradient of the sound- 
speed profile, a reasonably accurate solution is possible. There are some limitations 
involved with the use of ray acoustics and the model-based phase weights. These 
limitations will be discussed in the next section. 

C. LIMITATIONS OF RAY ACOUSTICS SOLUTION 
1. Turning Points 

A turning point is that position along a ray path propagating through an 
inhomogeneous medium at which the angle of propagation measured with respect to 
the positive Y axis, P(y), is equal to 90 degrees. At this point the origination of the ray 
path becomes ambiguous to a receiver using the localization technique described in this 



32 



thesis, because there is no way of knowing how many turning points the acoustic 
signal has passed through. The turning point will cause a transmitter that is below 
(above) the receiver to appear to be above (below) the receiver. Figure 2.5 illustrates 
these two possibilities. In the case of receiver one in Figure 2.5, a turning point has 
occurred between the transmitter's location and the receiver's location. The theory 
presented in this section would result in a calculated line-of-sight similar to that shown 
in Figure 2.5. The acoustic signal passes through two turning points prior to reaching 
receiver two, and the resulting line-of-sight calculation would indicate that the 
transmitter is at a depth below receiver two. 



Speed of Sound 

0 ► 




The turning point ambiguity problem is not necessarily very restrictive, 
depending on local sonar conditions. Table 1 lists the location of turning points in 
terms of AY and AR between the transmitter and receive array. The values in fable 1 
were calculated by assuming the values for (5(y 0 ), c(y 0 ), and g shown in Table 1, and 



then using equations 2.11 and 2.17 with P(y^) = 90°. These results show that for 
most angles of transmission the floor of the ocean would be reached prior to the signal 
reaching a turning point. Even at angles of transmission greater than 60 degrees the 
ranges to a turning point are quite large. 



TABLE 1 

DEPTH AND RANGE TO TURNING POINTS FOR A POSITIVE 

GRADIENT 



O * 


AY (km) 


AR (km) 


10° 


412.913 


492.043 


O 

o 


166.935 


238.343 


o 

O 


86.765 


150.276 


o 

O 


48.241 


103.424 


o 

O 


35.921 


86.765 


o 

o 

'O 


13.449 


50.063 


o 

o 


5.570 


31.582 


oo 

O 

o 


1.336 


15.271 


o 

OO 


0.331 


7.592 


c(y Q ) = 1475 m/sec 


g = 0.017 sec" 1 





If the transmitter and receive array are located in the negative gradient 
portion of the sound-speed profile as shown in Figure 2.5, the situation becomes much 
more restrictive. Here the transmitter must transmit in the upward direction to reach a 
turning point, as opposed to the downward transmission assumed in Table 1. Table 2 
contains the results of calculations for the turning points in this region. In this case, 
the angles were only varied from 91 degrees to 100.8 degrees in order to place the 
turning point within the negative portion of the sound-speed profile of Figure 2.5. 
Even with the higher magnitude gradient used in Table 2, ranges of several thousand 
meters are achievable prior to the turning point. Note that all distances in Table 2 are 
in meters, whereas those listed in Table 1 are in kilometers. 
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